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Abstract 

We study, on a d dimensional hypercubic lattice, a random walk which is homo- 
geneous except for one site. Instead of visiting this site, the walker hops over it with 
arbitrary rates. The probability distribution of this walk and the statistics associated 
with the hop-overs are found exactly. This analysis provides a simple approach to the 
problem of tagged diffusion, i.e., the movements of a tracer particle due to the dif- 
fusion of a vacancy. Applications to vacancy mediated disordering are given through 
two examples. 
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1 Introduction 



The venerable problem of a random walk on a regular d-dimensional hypercubic, infinite 
lattice continues to generate considerable interest from both novel quantities as- 

sociated with the simple walk and new variations of the walk itself. In the simplest case 
(which is called Polya walk M), the walker moves to one of the nearest neighbor sites 
with probability l/2d at each time step. One of the variations involves a "taboo" site, 
to which the walker may never visit. The question that is usually asked refers to the 
probability || , Q ( "taboo probability" ) of the walker visiting a certain site on the n-th 
step without having visited the taboo site on any of the erlier steps 1,2, ...,n — 1. An 
alternative way to think about the taboo site is to consider it as an irreversible trap, so 
that the above question translates into a question on the survival probability . A slightly 
different situation is when the walker arriving at a neighboring site of the taboo, is let to 
remain stationary, with probability l/2d, instead of moving into the forbidden site. All 
these variations belong to the chapter of lattice walks with "defective" sites in the theory 
of lattice walks [||. In this paper, we investigate a further case of walks with a defective 
site, namely, the walker hopping over this special site. Since it is no longer "taboo" (in the 
above sense), it will be referred to as the "hop-over" site. We write the Master equation 
for such a walk, with 2d arbitrary hop-over rates. Thanks to translational invariance, we 
may choose to locate the hop-over site at the origin. In Fourier space, we obtain a closed 
expression, in terms of the inverse of a matrix, for the probability distribution generating 
function, given a particular initial position so- 

Similar to the derivation of "taboo probabilities", the straightforward way to find 
the hop-over probabilities would rely on a study of return probabilities and first passage 
times. In this approach, probabilities of returning to the neighborhood of the origin 
without hopping over it are exploited in infinite sums over individual hop-over attempts. 
A similar method was used for solving the related problem of tagged diffusion || Our 
present approach is different, considerably simpler and more compact in terms of formulas. 
Instead of infinite sums, only the inverse of a certain 2d x 2d matrix needs to be computed. 
In d = 2, where this walk displays the most interesting characteristics, the inversion is 
quite simple. 

We are also able to keep track of both the direction and the frequency of the hop-over 
attempts. As a result, an interesting application is the diffusive behavior of a tagged 
particle g 0]. Such a tag may also be thought of as a "passive walker" , namely, a particle 
which remains inert until it is forced to exchange places with a normal, "active" random 
walker. In metals or alloys, an impurity atom can play the role of a tag, while a vacancy 
diffuses as a typical walker ||-|15]. We will discuss how to map the hop-over problem 
to the tagged walk. Finally, from the frequency of hop-over attempts, we can compute 
the distribution of "hits" received by the tag. We present two related applications: the 
first studies the rate of disordering of the surface of an A-B alloy, with a certain type 
of interactions, due to the wandering of a surface vacancy, while the second investigates 
the asymptotic time evolution of the disorder caused by a single Brownian vacancy in an 
initially completely phase seggregated, two-species system. 
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This paper is organized as follows. In the next section, the model and the Master 
equation are carefully defined (Sections 2.1 and 2.2) and the explicit solution is presented 
allowing for the extraction of the statistics on the hop-over events (Section 2.3). Section 3 
presents the relationship to the tagged diffusion problem with its complete solution through 
the techniques of Section 2., along with two associated applications presented in Sections 
3.2 and 3.3. The last Section is devoted to a summary and possible generalizations. 



2 Random Walk with a Hop- Over Site 



For completeness, we devote the first subsection to definitions, notations, and some well 
known properties of a random walk on a lattice. 

2.1 The simple random walk 

On a infinite (i-dimensional hypercubic lattice, the sites are labeled by s while the set of 2d 
lattice vectors is denoted by {a}. A walker, performing a pure random walk of Polya type 
(U, Q moves from a site to one of the 2d nearest neighbor sites at each time step (i.e., from 
s to s + a) with probability p. If 2dp < 1, the walker remains stationary with probability 
1 — 2dp. Given an initial position so we are interested in the probability distribution of 
finding the walker after n steps at site s : P^ (s|so) . Here, the superscript F reminds us 
that this is the distribution for "free" diffusion. The equation governing this walk is 



P n+l( S \ S 0) ~ P n ( S \ S 0) 

with the initial condition: 



p^[P r f(s + a|so)- J P„ F (s|s ) 

{a} 



P (s|s ) = <5 s , So 



(1) 



(2) 



Due to translational invariance, the solution can be obtained simply in terms of the gen- 
erating function in Fourier space: 



P(s|s ;£) = £P„(s|s )£ n , |£|<1, 

n=0 

P(k;£) = £P(s|s ;£y k - s ■ 

s 

Applying (||) to s = in (||), we obtain a useful identity 

&^P F (a|s ;0= [P F (0|s ;£)-<5 ,s ] + £(2dp - l)P F (0|s ; 

{a} 

Continuuing, we substitute (@) into (|) and (|) to arrive at 



(3) 
(4) 

(5) 



P F (k;0-e 



ikso 



P F (k; = -pJ2 [1 - cos(k-a)] P F (k; 

{a} 
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The solution is trivial: 

P F (k;0 = G(k,Oe rkso (6) 

where 

G(k,o = |i-e+£ P g[i-cos(k.a)]J (7) 

is the well known propagator for free diffusion. 
The inverses to (|j, ||) are given by 

P(s|s ;£)= / e- iks P(k;e) (8) 

and 

P n (s|so) = jfr^(s|so;0 > (9) 

where 

d d k , / 1 /d£ 

— 7T C 

and C is a suitable (counterclockwise) contour around £ = 0. Thus, the solution to the 
simple random walk can be written as 

P F (s\s ; = J e~ ik < s - s ^G(k, £) (10) 

and 

P F (s\s ) = J^^C n e- lk < s - So) G(k,0 . (11) 

In subsequent sections, several probability distributions will occur frequently. For 
convenience, we summarize their properties here. Define 

t = P F (0\0;0, (12) 

u = P f (8l\0;£) , (13) 

h = P F (a\- a -0, (14) 

v = P F (a\h; £) , b / ±a , (only for d > 2) (15) 

all of which are known functions of £. In particular, is just the generating function for 
the return probability of a pure random walk. The definitions (|l2|), ([l4| ) and (|l^) allow 
us to write 

P F (a|b; £) = v + (t - v)S aih + {h- v)S a ^ h , (16) 

where both a and b are lattice vectors (nearest neighbors of the origin). Applying (|5|) to 
the simplest random walk (p = l/2d), we find relations between t,u,h, and v ||, e.g., by 
choosing so = in (||) 

i = 1 + ^ , (17) 
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and (for d > 2) by choosing s = b 

£p[i + 2(d - l)u + ft] = u . (18) 

Finally, since we will be interested in late times (n — ► oo) corresponding to the limit 
£ — » 1~, we note the following well known asymptotic behavior of t Q: 

f 1/^2(1-0 ford=l 

^ ^ ln (i^) ford=2 (19) 

[ const. for d > 2 

From dl7|) , we conclude that u behaves the same way. 



2.2 Walks with a hop-over site 

Without loss of generality, let us place the hop-over site at the origin. The probability 
for the walker to hop from —a to a will be denoted by p a . By keeping this rate different 
from p, we will be able to log the different hop-over attempts. Clearly, P r ^(s|so), the 
distribution with such a site, will be a polynomial in the various p a 's. Meanwhile, the 
co-efficient of p v a will be associated with the subset of those walks which have hopped 
over the origin (from —a to a) f a times. We will return to these considerations in more 
detail below. 

Again, we let the walker start at so, which is not the origin, i.e, 

P(f (s|s ) = 5 s , So , and s ± . (20) 
The subsequent evolution of P H is governed by the following Master equation: 
P? +1 (s\s ) ~ Pn(s\s ) = P (l - 5 a , ) ]T [P* (s + a|s ) - i^(s|so)' 

{a} 

+ E 6 *,» [Pa^(-a|s ) - (p_ a - p) P% (a|s )] (21) 

{a} 

the first term on the right shows that the walker never visits the forbidden site (0) and 
that, away from the neighborhood of the origin, it performs a simple Polya walk (a "free 
walk"). The latter term describes the possibility of hop-over, when the walker finds itself 
on a nearest neighbor site of 0. (See Fig.l.) 

Going over to (k, £ ) space by (f|J|) , this equation becomes 

G- 1 (k;0P H (^0-e ikso = -^P H (a|s ;£) 

{a} 

+£E {p a P H (-*\so;0 - (p- a -p)P H (a\s ;0] e ik ' a • 

{a} 
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This representation clearly displays the effects of the "defect" (associated with the hop- 
over site), since the left hand side consists of the terms for the pure random walk, only. A 
more elegant form would be 

G-^k; £)P*(k; - e tks ° =£E r ( k < a)P H (a|s ; (22) 

{a} 
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Figure 1: Random walk with a single hop-over site. The diamond symbolizes the hop-over 
site located in the origin and the circle represents the vacancy. 



where 



T(k, a) = p (e ika - l) + p_ a (V ika - e ika ) (23) 
Now, the right hand side of (e3) can be regarded as an extra inhomogeneity for the solution: 



P"(k;0=G(k;0 



(24) 



e !ks "+^r(k,a)P H (a|s ;0 

{a} 

To find the solution explicitly, we must determine the 2d quantities P^(a|so; £). This 
can be done by exploiting (]§) and (| 



i^(a|s ;0 



e- lk - a P^(k; 



-ika 



G(k;0 



°+eE r ( k ' b ) pH ( b l s o;6 



{b} 



P F (a|s ;0+E^ / e- ik - a G(k;Or(k,b)P^(b|s ;0 (25) 
{b} Jk 



Linear in P^(a|so;£), equation (|25| ) can be solved: 

P^(a|s ;0 = E L a,b(0^ F (b|s ;0 



(26) 



{b} 
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where L is the inverse of the 2d x 2d matrix with elements: 



a,b — °a,h 



-ska 



G(k;Or(k,b) 



(27) 



Once P H (&\so;£) are known, they can be substituted back into (24) and the explicit 
solution can be obtained: 

P H (k; = P F (k; + E £ G ( k 5 Or(k, a)L aib (e)P F (b|s ; • (28) 

{a,b} 

Using we may return to the lattice: 

P H (s\s ;0 = P F (s\s ^)+ W£/ e- lk - s G(k;e)r(k,a) 

{a,b} L Jk 



L a , b (OP F (b| So ;0 • (29) 

Making use of Eq. ( |23| ) the elements of L _1 become expressed solely in free walk terms: 
(L-%, h = 6 a , h - i { P - h P F (a\ - b; + (P - P-b)P F (a|b; £) - pP F (a\0; £)} . (30) 
The integral in the square bracket of Eq. ( pS[ ) is computed similarly: 

£ / e- iks G(k;e)r(k,a) =^- a P F (s| -a;e)+C(p-p- a )P F (s|a;e)-ep J P F (s|0;e) (31) 



Equations ( |29[]3l| ) express the generating function of the hop-over walk in terms of the 
hop-over rates {pa} and the generating function characterizing the free walk. It may be 



worthwhile to interpret the second term in (psj) physically. From (27) L can be expanded 
as a power series in T, which accounts for the difference between a free passage through 
the origin and an hop-over. Thus, the n th term in this T series in ( |29|) can be thought of 
as paths which have n encounters with the hop-over site. This connection will be explored 
further in the next subsection. 

Here, we end with a closer look into the case with isotropic hop-over rates, i.e., 

P a =p' , 

but with p' not necessarily being equal to p. Then, eqns. (26-31) reduce to 

P H (s\s ;Z) = P F (s\s -0+t £ [p>P F (s\- a ;Z) + (p-p')P F (s\ a -0 

{a,b} 

-pP F (s|0;£)lL a , b (£)P F (b|so;0 



while our task is to invert, 

(L _1 ) a ,b = <5 a , b - £ \p'P F (a\ - b; £) + (p - p')P F ( a \h; £) - pP F ( a \0; £) 



(32) 



(33) 
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where a and b are nearest neighbors of the origin. Due to isotropy and homogeneity, this 
matrix reduces considerably. The last term in the [...] brackets above is seen to be a 
constant (^), — pu, for all matrix elements. Using (|l6|), we write (^) as 

(L^Kb = Zp(u -v)+ 5 a , b [1 - i (p'(h - t) + pit - v))] 

-5 a ^U[p'(t-h)+p(h-v)] (34) 

Since this matrix is of the form 



A + BS a>h + CS a - h , 
it is easy to check that its inverse is also of this form, namely, 



J a,b 



A' + B'5 a u + C% 



a,-b 



1 



B 2 -C 2 



( A(C-B) 
\2dA + B + C 



+ B6 a , h - CS a _ b 



(35) 

(36) 
(37) 



With the help of (|34j37| ), we can carry out the matrix multiplication in (p2\). Note that a 
quite a few sums over {a, b} decouple. For example, using (|5|) and P F (s|s'; £) = P F (s'|s; ^) 
, we have 



£ J2 {p'P F (s\ - a; + (P - p')P F (s\a; ~ pP F '(s|0; 

{a} 



(l-OJ^(s|0;0 



The explicit expression for our walker to move from so to s, with isotropic hop-over rate 
p', is 

P H (s\s ;0 = P F (s\s ; e) - [B' + C" - A' (1 - £) /$p] [1 - $ + 2^p] P F (s|0; 0^ F (0|s ; 
H [B'p + (C - B>) p'] P F (s\*; Oi^(a|so; 

{a} 

H [C'p + (-B' - C')p'] P F (*\ ~ a; 0^ F (a|s ; (38) 

{a} 



where A', P/, and C" are to be read off from (34-37). Except for the first term (which 
stands for the 'free' case), this expression shows the non-trivial effect of a single hop-over 
site on the random walk explicitly . The various terms can be easily interpreted: the 
second term removes the (free) walks which pass through the origin, while the last two 
takes into account walks which land on a neighbor. If we take the usual continuum limit 
(the lattice spacing a — ► 0, unit time-step r — > 0, such that a 2 /r = 2d) of this expression, 
the effect becomes vanishingly small. This result is perhaps not surprising, since a single 
site cannot affect the properties of the random walker in the large distance, long time 
limit. However, buried in this approach is a non-trivial question, namely, the statistics 
associated with the hop-overs (how often and in which direction does the walker hops over 
the origin). That is the subject of the next section. 
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2.3 Statistics of hop-overs 

Next, we turn to a study of the statistics of the hop-overs. Note that, due to the longer 
jumps (across the origin), the walker suffers an extra displacement when compared to the 
simple free walk. We define the variable p as the negative of this extra displacement, for 
reasons that will become obvious later. We will also be interested in the number of times a 
hop-over jump has occurred. Thus, we define <^(s, p|so) as the probability that, starting 
from so, the walker arrives at site s after n steps, 

• having performed v hops over the origin and 

• suffering a total "extra" displacement equal to —p. 

As usual, it is more convenient to work with generating functions in Fourier space. 
Accordingly, we trade the variables (s, p ) and (n, u) for their conjugates: (k, k) and 
(£)0- We can write: 

oo 

*(k,K;6C)= E Erre i(k - S+K - p) C(s,p|so) ■ (39) 

n,u=0 s,p 

Once $(k, k; £, Q is known, the inverse transforms will lead us to the distribution itself: 
<(s,p|s )=/ / r^e- i(k ' s+K - p) *(k,K;£,C) ■ (40) 

To arrive at an expression for ^(k, k; £, £), we first study a more detailed distribution: 
0^ a ^(s|so), i.e., the probability of the walker makes exactly v a hops from —a to a (for 
each a). Here {f a } denotes the set of 2d numbers associated with the different directions 
of hop-over. Since each hop-over may occur only with probability p a , and the events are 
independent, we conclude that ^^(slso) must contain the factor Il{a}Pa a > where II{ a } 
is product over all the 2d nearest neighbors. 

Meanwhile, the extra displacement is just —p = J2{a} azy a> so that 

C(s, p|s ) = E 4 M (s|so) <y„ >E{a} „ a 5„ p>E{a} aUa . (41) 
{"*} 

In terms of generating functions, this expression becomes: 

<l(k,K;£,C) = E * K} (k;0nC* a e-^ a (42) 

{a} 

where $^ a ^(k;^) is the generating function and Fourier transform of ^l^^ (s\so) . 

Armed with these considerations, we interpret the expression ( f28| ) physically. Expand- 
ing L in a power series in the p a s, we see that P H itself is a power series in these rates. 



9 



Naively, it is tempting to identify these co-efficients with 0^ a ^(s). However, there is an 
implicit probability for the walker to remain at the site —a : 

(l-p a -(2d-l)p) , 

so that the co-efficient of p a also includes paths that avoid the origin. Another way to 
see this difficulty is through the original Master equation (21), in which a particular p a 
appears twice, once in the "gain" of P H (a) and once in the "loss" of P H (— a). Of course, 
both of these terms represent the same jump. Therefore, to avoid double counting, we 
can label one of them by q a and, only at the end, set q a = p a . Let us emphasize that, 
with q a 7^ p a , total probability will not be conserved. But this is precisely the trick for 
distinguishing making a hop-over jump from not making one. 

For definiteness, let us choose to label the "gain" term by q a . A little care leads us to 
a modified version of ( 23 ) : 

r(k,a) = p(e jka -l) +(? _ ae - lka -p_ ae ika 

= q _ ae -^ + (p-p__ a ) e ^-p (43) 

from which we can find 

f H (k;0=P F (k;0+ E ^(k;Or(k,a)I a>b (e)P F (b;0 . (44) 

{a,b} 

Note that the overline over a quantity symbolizes the dependence on both sets of rates: 
{q a } and {p a }. 

With this trick, we can associate the co-efficient of (q a ) u * , in a power series expansion 
of P H , with walks that include v a hop-overs (jumps from —a to a). Now, we must pick 
out these co-efficients and then multiply them with the correct weight: (pa)"*- This is 
accomplished by applying p v a & § q a v& for each a, so that 

*^>(k;e) = Il/ (Pa/fe)" a ^ H (k;0 • (45) 



{a} 



Substituting Eq. ( [451) back into (|42|), we see that the dependence on the variables v a 
factorizes. Performing the sums over each v a , we arrive at: 

<&(k,K;£,C) = II^/ %3i^(k;0 • (46) 

f\ 2m J q a -p a (e lKa 



i.e. 



$(k,K;£,C)=P H (k;0 



(47) 

<?a=PaC e 



In other words, <J>(k,K;£, () is nothing but P H (k;( i ) with all the q a s simply replaced by 
p a (e~ lKa . 
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Prom expression for P , we write explicitly: 
$(k, k; e, C) = ^ F (k; + E ? G ( k ; r (k, «, a; C)La,b(^; £, C)^ F (b|s ; , (48) 

{a,b} 

where 

r(k,K,a;C) = P- a 

and 

(IT 1 (k; £, C))a,b = <5 a ,b - e {^bC^ b ^ F (a| - b; f) + (p - P - h )P F (a\b; £) - pu) (50) 

Due to the extra factor £e* K ' b , it is difficult, even for the isotropic hop-over walk, to 
invert L in arbitrary d. However, since most of the interesting effects occur in d = 2 (where 
the return probability approaches unity for large times), the 4x4 matrix can always be 
inverted. 

2.4 A hidden pure random walk 

Finally, if we are interested in the simple hop-over problem in which all the rates are 
identical, then we simply set all p a 's to p. In this case, we should also retrieve the pure 
random walk if we subtract the extra displacements due to the hop-overs. In other words, 
let us consider the probability of finding the walker at s— (— p), regardless of how many 
hop-overs occured or their directions. Thus, we study 

oo 

*n(r|so) = EE 'Wp E ^n( s > Pl s o) • 

s p u=0 

Again, it is easier to study the generating function in Fourier space: Using k as the variable 
conjugate to r, we see that the above sums result in: 

P(k;0 = E^'T^nMso) 

n,r 

= $(k,« = k;£,C = l) • 

Setting in (|49| ) p_ a = p for each a, k = k and ( = 1 we see that T(k, k = k, a; £ = 1) = 0. 
Thus, from (f44|), we obtain 

P(k;£)=P F (k;0 , (51) 

justifying our expectation that, if we subtract the displacements due to hop-overs, the 
simple random walk re-emerges. 



-i(k— re)-a 



tka 



+ p[e 



,ika 



(49) 
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3 Relation to Tagged Diffusion and Some Applications 



The problem of tagged diffusion has been solved previously (in two-dimensions), ||] (in 
(^-dimensions) , using a "direct" approach, i.e., by keeping track of all the possible ways 
for the tag to be displaced by a random walker. In this section, we will first show that the 
above analysis can be applied to provide a new approach to this venerable problem. As 
examples of applications of our results, we devote the latter subsections to the dynamics 
of two particular physical systems. 



3.1 The passive walk or tagged diffusion 

For completeness, let us describe tagged diffusion in terms of a "passive" random walker. 
Returning to our infinite d-dimensional hypercubic lattice, we place two walkers, one 
active and one passive. Respectively, we label them B (for Brownian) and T (for tagged, 
or tracer, particle). The former performs the simplest random walk, as in Section 1.1 with 
p = l/2d. When B attempts to move to the site occupied by T, the two simply exchange 
places. Thus, T does not move except when it is "kicked" by B. In this sense, the motion 
of T could be called "passive" and will be referred to as a "Brownian driven walk" . The 
mathematical properties of this pair of walkers are contained in the joint probability 

C(r,p|ro,0) (52) 

for finding B at site r and T at site p, on the n-th step of B and the j/-th step of T. Note 
that, v just represents the number of "kicks" T received from B. The last arguments refer 
to the initial condition, i.e., T being at the origin and B at site ro ^ 0. 

From the description above, it is clear that, relative to T, B is performing a random 
walk with a hop-over site. Thus, 

C(r,p|r ,0) = C(s,p|so) (53) 

of the previous section, provided ro = so, and 

s = r-p (54) 

which is just the position of B relative to T. Of course, we expect the results following this 
approach to be identical to those from the more "direct" approaches in || and . To be 
brief, here we will present only a certain projection of the distribution <J>^(r, p|ro, 0), and 
compare the results with those in Ref. ||. 

Let us focus on the probability that T received v 'kicks' during n -steps of B, regardless 
of the final locations of the two walkers. Denoting this quantity by <^(ro), it is: 

<(ro) = ^<(r,p|ro,0) = ^^(s,p|s ) . (55) 

r,p s,p 

From (39) and (f47|), we see that the generating function 

o = E err 



12 



is given by 



C) = *(0,0;& 0=^(0,0; &C) . 



(56) 



Since B performs a random walk, we set p a = p everywhere. Thus, from eqn. (|4q), we 
have: 

1 



l+e E r(0,0,a;C)L aib (0;C,C)P F (b|s ;0 . 

{a,b} 



1-C 

Using (H), (|§ and (|l|), we find: 

r(0,0,a;C) = -p(l-C) 

and 

(ir 1 (o; e, c))a, b = ep (« - Cv) + [i - ^cp(^ - «)] <y a ,b - ecp(* - «)*a,-b 

Following we write L a b in the form 

L a ,b = A' + B'^ a ,b + C"^a,-b , 

so that 



1 - £p(1 - C) (2dA' + 5' + C") E ^ F (b|s ; . 

{b} 



(57) 

(58) 
(59) 

(60) 
(61) 



Now, to simplify this some more, we make use of (|5|) and carry out the algebra for A', A, 
etc.: 



1 



(i-0 



1 — e L 1 + £u + CCP ((2 - 2d)v -h-t)) 



P F (0|s ;0 



Finally, exploiting (17, lq), we write a compact form for the generating function: 

1 



I - ai : C) . P F (0|s ;g) 



(62) 



(63) 



# - c* + c 

Before inverting this result to obtain let us explore the scaling properties of £), 
in the limit £, £ — > 1~. For d < 2, t diverges according to (|l9|), but (1 — £)t — ► 0. So 
does P F (0|so;^) for any so, with the result that P F (0\so; /t — ► 1. Keeping the leading 
non-trivial orders, a simple scaling function emerges 



(1-60(^0 



i + x' 



where 



x= (l-C)t 



(64) 



(65) 



is the scaling variable. If we trace the origins of (|19|), we can use an explicit and more 
general form: x = (1 - <") (1 - £)( d - 2 )/ 2 . 

Returning to the inverse transforms, the one with respect to £ is trivial. Keeping in 
mind (17), we find that (~ u 4>(£,, £) is identical to eqns. (34) and (35) of Ref. || obtained 
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via a completely different method. Since t is a non-trivial function of £, the transform 
with respect to £ cannot be carried out explicitly. Thus, we leave the final result in the 
form of an inverse transform: 
v > 1 : 



n- 



(l + tt-Wl-l/t)"- 1 P*-(0\r ;0 , 



u = 0: 



#(ro) 



t-P F (0|r ;0 



(66) 



(67) 



These are the explicit formulae, in any d, for the probability that the tag particle (passive 
walker) moves v steps while the vacancy (active walker) takes n steps, starting atro from 
the tag. Of course, for v,n —* oo, the scaling analysis above implies that 



v ~ y/n and 



'Inn 



in d = 1 and 2, respectively. Above two dimensions, the return probability of the random 
walker remains less than unity as n — > oo, so that <\f n decays exponentially in v. 



3.2 Vacancy mediated disordering of an A-B alloy with extreme anisotropy 

In some binary alloys, an ordered state consists of A and B atoms occupying alternate sites 
on a cubic lattice. Under appropriate conditions, an atom cannot move unless a vacancy 
comes in contact and exchanges places with it (also coined as the 'vacancy-mechanism') 
f8|- fl5f| . Therefore, each atom can be regarded as a passive walker, while the vacancy plays 
the role of an active walker. More precisely, consider a monolayer, composed of A and 
B atoms in a square lattice, adsorbed on a substrate which interacts strongly with these 
atoms in such a way that the ground state (of the monolayer) is a simple checkerboard 
configuration (antiferromagnetic state, in the spin language). To simplify the problem 
further, we suppose that the intral&yer interactions are negligible. A similar arrangement 
may be realised by a binary alloy in a NaCl structure on a simple cubic lattice, with 
extreme anisotropic interactions. If the interactions are much stronger along the c-axis 
and we are focusing on a (0,0,1) surface, the surface atoms will experience a much larger 
interaction with the bulk than with their neighbors on the surface. For clarity, we will 
distinguish the surface atoms from bulk ones, through labeling the former by A/B and 
the latter by a/b. The ground state, is shown in Fig. 2a, with only Ab and Ba bonds 
(and, of course, AB and ab ones). Given our assumptions, Bb or Aa bonds, considered 
as "mismatches", will be more costly. Thus the movement of the vacancy can cause such 
excitations (Fig. 2b). In a typical system, the vacancy will not perform a pure random 
walk, since its movements will be governed by these excitation (or de-excitation) energies. 
However, if subjected to sufficiently high temperatures, we could neglect this dependence 
and approximate the vacancy by a Brownian particle. 

Focusing on such a vacancy wandering in our monolayer, we can investigate the evo- 
lution of the excitation energy through the creation of the mismatches. In particular, the 
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a) 




\ 

'mismatch' 



Figure 2: A binary alloy model with extreme anisotropy. The configuration a) "A on top 
of b" or "B on top of a" is energetically favorable to, e.g., the configuration b) "B on top 
of b" which is considered as a mismatch. 



analysis above can be applied to answer the following question: after a long time (n>l) 
what is the expected excitation energy caused by this vacancy? Being proportional to 
the expected number of "mismatches", this energy shift can be found readily. Note that, 
though the discussion above is clearly based on a d = 2 surface adsorbed on a d = 3 bulk, 
our considerations will be applicable in any d. 

First, let the vacancy be located initially at the origin. Then we tag each atom simply 
by its initial location. Alternatively, using flBlj), we may use so (instead of — so) to label the 
particle uniquely. After the vacancy has taken n steps, atom 'so' has suffered v displace- 
ments with probability 0^(so), which is found in (|55|). If v is odd, there is a "mismatch" 
for that particle. The expected total number of mismatches is therefore expressed by: 

<^> = E' Eo^-Hi*)- (68) 

so u=0 

where the prime on the summation symbol means that so = is excluded. Such a sum 
for the quantities in eqns. ( p6| ) and ( |67| ) is trivial: 

P(0|s ;£) = £P(0|s ;0 - P(0|0;0 = - t . (69) 

So S ^ 

After summing over u, we find the desired result: 

i i - n r(e-i)t+i 



(Yn) 



e(i-fl' 



(70) 



To proceed, let us consider the large n limit, corresponding to £ — > 1 . Since 
(£ — 1) t — * even for cases when t diverges, the [.] bracket can be replaced by 1/ {2t — 1). 
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The n — > oo behavior can then be extracted by exploiting the Discrete Tauberian Theorem 
(See, e.g., Ref. @, pp. 118.). 

In d = 1, where we have a linear chain of atoms (adsorbed on a two dimensional bulk), 
we find that the expected number of mismatches grows like: 

(Y n ) ~ v/m/? . (71) 

V 7T 

It is interesting to compare this result to the expected number of distinct sites visited in 
an n-step walk, which is a well known quantity pL In d = 1, it grows like: 

(S n ) ~ \/n\ —, as n — > oo . (72) 

V 7T 



Thus, we find that 

(Yn) 1 



(73) 



(S n ) 2 

as n -> oo, i.e., half of the visited sites will have mismatches. These results are hardly 
surprising, since, at any time, only the string of atoms between the walker and its initial 
position are "mismatched". 

For d = 2, the physical case, the results are more interesting. With the help of the 
same theorem, we find: 

( ^~2hTM' as n "°°' (74) 

whereas [[jj 

( Sn ^ ~ w^T ' as oo , (75) 

in (8nj 

so that (^) proves to be valid in this case also. 

Though bulk materials in 4 or higher dimensions are not physical, we may nevertheless 
consider the number of expected mismatches in d > 3. Here, t (1) is finite, leading to 

(Y n ) ~ n —4 • (76) 

x 1 2t(l) - 1 v ; 



Comparing it with 



t(\ 

we find, instead of ([73|), a more "interesting" result: 

(Yn) tJX) _ 1 

(5 n ) ~ 2£(1) - 1 1 + i2(0) 



(Sn) ~ n , (77) 



as n — > oo (78) 



where i?(0) is the probability that the walker ever returns to its starting site, another well 
known quantity in the theory of random walks. Actually, eqn. ( |7"8| ) is in fact valid for any 
d > 1. For d < 2, the walk is recurrent, i.e., -R(O) = 1, reducing (fftl) to 1/2. Only for 
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d > 2, is the walks transient, where R(0) < 1. An alternative way to display ( |78"D is to 
quote the ratio of particles performing an even number of exchanges with the hole to those 
"visited" an odd number of times. Clearly, this ratio is simply R(0). These remarks give 
an intuitive explanation for (|78|). Since the random walk is recurrent in d = 1,2, every 
particle is sure to be revisited. Therefore, there should be as many particles visited an 
odd number of times as those visited an even number of times, on the average. However, 
for d > 3, a certain amount of the particles will be visited only once. Thus, the particles 
visited an odd number of times should be larger. 



3.3 Propagation of interfacial disorder 

As another example of how our results can be applied in disordering dynamics, we in- 
vestigate the asymptotic behavior of the disorder induced by a single Brownian vacancy 



1 16]. Since the most interesting case is d = 2, we will limit our study here. The initial 
configuration is a completely phase segregated system, i.e., an infinite square lattice, with 
the upper half plane filled with one type of particles (white) and the lower half with the 
other type (black). A single vacancy (the active random walker) is placed in the white 
region at the interface between the two half planes (Fig. 3), and labeled as the origin of 
our coordinate system. As the vacancy wanders, particles will be drawn into the opposite 
phases, leading to disordering. In the following, we will show that, after n steps taken by 
the wanderer, the "disorder-profile" parallel and perpendicular to the interface scales as 
^fn and y/lnn , respectively. The scaling functions are combinations of exponentials and 
a modified Bessel function. 

Let us focus on the black particles. One definition of disorder is the density of black 
particles in the upper half plane. Thus, we seek 

the probability of finding a black particle at location s, after n steps of the vacancy. Since 
each particle performs a "passive" random walk independent of the others, we may simply 
track the movements of each black particle and sum over all those which reach s. As in 
previous sections each black particle can be labeled uniquely by its initial position so- 

Since we are not interested in neither the location of the vacancy nor the frequency of 
a particle being kicked it is clear that we should sum over the vacancy's position and the 
frequency of the hits the black particle receives. Summing over all the black particles, we 
see that 

(-) (-) 
K(*) = E <M S " so | -so) = E E E <( r > s - s ° I "So), (79) 

so so r u=0 

where denotes summation over the so's in the lower half plane only , and r is the 

walker's position after n steps relative to the site so- In terms of the generating function 
(c.f. eqns. ( |39l) and (p0|)), we have 

$>;£)= / ° )K *(0,K|-8o;e,l) . (80) 



so 
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Figure 3: Initial configuration of a sharp interface between two different species (transpar- 
ent and filled squares). The vacancy (large empty circle) is located initially at the origin 
of the coordinate system. 

where the sums over v and r have been carried out and the explicit dependence on so 
restored. Referring to eqns. (f48|)-([5Ci|), we have 

*(0, k\ - s ; e, 1) = 7^7 | 1 + £ E f (°> K > a ; 1 )^a,b(^; S, l)P F (h\ - s ; 1 , (81) 

1 " * I {a,b } J 

r(0,^a;l)=p(e iK - a -l) (82) 

and 

(L- 1 (K;Z,l)) a , h = 5 a , h -&{e iK - h P F ( a \-b;0-u} . (83) 

At this point, we will restrict our attention to d = 2 only. Clearly, there is no interesting 
behavior in d = 1. On the other hand, the random walk being transient in d > 2, disorder 
will be confined in the large time limit. 

In two dimensions, the matrix (^) is 4 x 4 and must be inverted laboriously. The cal- 
culations involved from (^0|) to ( |83| ) are reasonably straightforward, but extremely lengthy. 
We shall present only the final result obtained in the n> 1 limit. 

For fixed s = (si, S2), we obtain 

* , x 1 1 + 2«2 / ; 7 1 hi In n „„ / 1 \ , 
f (s) = — - Jtt(tt - 1 — = - + O . 84) 

It is clear that the complete disorder is the final state, given by (ft* = 1/2, i.e., "complete 
gray". As expected, this value is approached from below if S2 > 0, and from above if 
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S2 < 0. Since the frequency of "kicks" scales as \/lnn, we see that the decay follows 
1/Vlnn rather than the typical decay of a random walk, i.e., 1/y/n. 

Instead of fixed s, we seek a scaled distribution. For this we note first that disorder 
along the interface should arise relatively quickly, since it depends only on the presence of 
the random walker, which wanders afar as y/n. On the other hand, disorder in the vertical 
direction (from the origin) relies entirely on the wandering of the passive walkers, so that 
it will occur at the Vlnn time scale. The appropriate scaling variables turn out to be 

x = 2s\j\fn and y = 2.S2 \Jtt(tt — 1)/ Inn 
while the final result is 

W n {x,y) - Q(-y)} ln8n = sgn(y)K (\x\)e-^ + O (^=) (85) 

Vvmn/ 

where sgn{y) denotes the sign of y and Kq is the modified Bessel function. Note that 
the square bracket on the left hand side is a measure of the disorder, since © is the 
step function that represents the initial probability distribution for finding black particles. 
Finally, since Kq{z) — > ^ir/2z exp (—z) for large z, we see that the decay in both directions 
are dominated by exponentials (in the scaling variables x and y). 

4 Summary and Outlook 

We have presented a study of a class of random walks on a d dimensional hypercubic 
lattice in which the walker hops from site to nearest neighbor site with one exception. It 
hops over a particular site. In case the hop-over rates are isotropic, the full probability 
distribution was found explicitly. We have also investigated the statistics associated with 
the hop-overs, so that we can find both how often and in which direction the hop-overs 
occur. The latter study can be readily applied to the behavior of a tagged particle, which 
moves only if a mobile vacancy were to exchange places with it. All previously known 
results of tagged diffusion are easily recovered in this simpler, novel approach. Finally, we 
showed how this study can be applied to two physical examples of the vacancy mediated 
disordering process. 

Though we have focussed only on infinite systems, our methods are readily general- 
izable to finite (periodic) lattices. In three or higher dimensions, such a generalization is 
crucial, since the random walk is transient, so that an infinitesimal density of vacancies 
can not give rise to system-wide dynamics. One way to estimate the effects of finite den- 
sity is to consider finite systems with periodic boundary conditions. Another interesting 
generalization is the biased random walk. With a bias, there would be little reason to 
study hop-over walks on an infinite lattice, since the walker never returns to the hop-over 
site. By contrast, on a finite lattice, the walker will "run into" the special site periodically. 
In the steady state, i.e., for times large compared to the traverse time, we can expect an 
asymmetric distribution around the "defect" . If the hop-over rates are also biased, say, in 
the opposite direction, it may be possible for an effective binding to occur. This problem 
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can be mapped into a limiting version of the biased diffusion of two species introduced 
sometime ago p7| . It would be very interesting to examine the walker-defect distributions 
and check if long range correlations [18] also appear. 
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